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ABSTRACT 

We present the results of the first multiwavelength observing campaign on the high-mass X-ray 
binary LS I +61303 comprising observations at the TeV regime with the MAGIC telescope, along 
with X-ray observations with Chandra, and radio interferometric observations with the MERLIN, 
EVN and VLBA arrays, in October and November 2006. From our MERLIN observations, we can 
exclude the existence of large scale (~ 100 mas) persistent radio-jets. Our 5.0 GHz VLBA observations 
display morphological similarities to previous 8.4 GHz VLBA observations carried out at the same 
orbital phase, suggesting a high level of periodicity and stability of the processes behind the radio 
emission. This makes it unlikely that variability of the radio emission is due to the interaction of an 
outflow with variable wind clumps. If the radio emission is produced by a milliarcsecond scale jet, it 
should also show a stable, periodic behavior. It is then difficult to reconcile the absence of a large 
scale jet (~100 mas) in our observations with the evidence of a persistent relativistic jet reported 
previously. We find a possible hint of temporal correlation between the X-ray and TeV emissions and 
evidence for radio/TeV non-correlation, which points to the existence of one population of particles 
producing the radio emission and a different one producing the X-ray and TeV emissions. Finally, we 
present a quasi-simultaneous energy spectrum including radio, X-ray and TeV bands. 
Subject headings: gamma rays: observations, X-rays: binaries, X-rays: individual (LS I +61 303) 
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1. INTRODUCTION 

LS I +61 303 is a high-mass X-ray binary consist- 
ing of a low-mass [M ~ (1 — 4) M Q ] compact ob- 
ject orbiting around an earl y type BO Ve star along 
an eccentric (e = 0.7) orbit (|Casares et al.l 120051 and 
refer ences therein). The modulation of both , the ra- 
dio (IGregorv fc Tavlor I [l978t IGregor yljgOOl and X- 
ray (jTavlor et al. 1 119961 : iParedes et al.lll997h emissions, 
display a period of -Pqj-1-, = 26.496 d, attributed to 
the orbital motion. The position of the maximum of 
the radio emission along the orbit, as well as its in- 
tensity, are modulated with a supero rbital period of 
P sup = 1667 ± 8 days (| Gregory! [2001) . LS I +61303 
is positionally coincid ent with an EGRET 7-ray source 
(Knif fenet al. I fl997h . Moreover, variable emission at 
TeV energi es has been recently d etected with the MAGIC 
telescope (| Albert et al. I l20"06aD . These authors found 
that the peak flux at TeV energies occurs at orbital 
phase orD ~ 0.65, while no high-energy emission 
is detected around periastron passage (4> or u ~ 0.23). 
From ~50 mas resolution radio image s of LS I +61303 
obtained with MERLIN, iMassi et all (|2004h suggested 
the existence of a precessing relativistic ((3 « 0.6) 
jet up to angular scales of ^0.1 arcsec, which led 
them to interpret LS I +61 303 within the framework 
of the microquasar scenario (jBosch- Ramon et al.1 120061 
iRomero et al. Il2007f). However, re cent VLBA imaging 
obtained bv iDhawan et al. I (|2006l ) over a full orbit of 
LS I +61303 has shown the radio emission to come 
from angular scales smaller than about 7 mas (projected 
size 14 AU at an assumed distance of 2 kpc). This ra- 
dio emission appeared comet ary- like, interpreted to be 
pointing away from the high mass star within a par- 
ticular scenario, and no relativistic motion, nor halos, 
nor larger-scale structures were detecte d at any phase 
of the orbit. Based on these findings, IDhawan et al. I 
( 2006) concluded that the radio and TeV emissions from 
LS I +61 303 are originated by the interaction of the 
win d of a young pulsar with t hat of the st ellar compan- 
ion (jMaraschi fc Treveslll98lUDubusll2006f K 

In this article, we present and discuss the results of a 
multi-wavelength campaign including radio, X-ray, and 
TeV 7-ray observations on LS I +61 303 , aimed at shed- 
ding light on the physical processes going on in the sys- 
tem, as well as yielding useful input for a detailed, time- 
dependent modeling of this relevant system. For this, 
we study the correlations between the different detected 
emissions in terms of their morphological, temporal and 
spectral features, both in the intra-day and day-to-day 
timescales. 

2. MULTIWAVELENGTH OBSERVATIONS 

Table Q] and Figure Q] summarize our observations of 
LS I +61 303 , carried out during October and November 
2006. In particular, we set up a simultaneous, multi- 
wavelength campaign on LS I +61 303 for the 25 and 
26 October 2006 using the MERLIN, EVN, and VLBA 
interferometers (radio), and the TNG (infra-red), Chan- 
dra (X-rays), and MAGIC TeV 7-rays) telescopes. Un- 
fortunately, bad weather conditions did not allow us to 
get any infra-red data, nor did we get useful data from 
the MAGIC telescope coincident in time with simultane- 
ously scheduled Chandra observations. The period 16-20 



TABLE 1 

Observing log for the multi-wavelength campaign 



Telescope 


Date 


UT 


range 


^orr/ 


Average flux 


MFRTJN 


26 Oct 


22:32 


— 04:30 


0.66 


34 84-1-0 93 




27 Oct 


10:30 


— 23:59 


0.68 






28 Oct 


00:00 


— 10:00 


0.70 


on fiQ-j-o 1 K 




Ifi Nov 


14:43 


- 23:59 


0.44 


71 41 _|_n 1 r: 




17 Nov 


00:00 


— 23:59 


0.47 


(U.I O^U. i- 1 *- 




18 Nov 


00:00 


- 23:59 


0.50 


78.98±0.21 




19 Nov 


00:00 


- 13:25 


0.53 


69.05±0.15 




20 Nov 


00:00 


- 11:24 


0.57 


47.14±0.25 


EVN 


26 Oct 


22:30 


- 04:30 


0.65 


33.29±0.27 


VLBA 


25 Oct 


21:30 


- 02:50 


0.61 


39.49±0.18 




26 Oct 


21:30 


- 02:50 


0.65 


32.69±0.14 


Chandra 


25 Oct 


22:13 


- 04:26 


0.61 


1.87±0.16 


MAGIC 


27 Oct 


00:50 


- 04:29 


0.66 


1.7±0.4 




16 Nov 


21:33 


- 01:04 


0.45 


<1.19 




17 Nov 


20:56 


- 01:00 


0.48 


<0.80 




18 Nov 


21:00 


- 22:19 


0.52 


<1.51 




19 Nov 


21:00 


- 22:00 


0.56 


<1.31 



a The orbital and super-orbital phases are com puted using 
MJD = 43366.275, P orb = 26.4960 flGreeorvll2002D . The super- 
orbital phase is </>sup = 0.4 for all the observation period. 
b The average measured fluxes are in mjy for radio observations 
(MERLIN, EVN, VLBA), 10~ n ergs cm" 2 s" 1 for the Chan- 
dra X-ray observations (absorbed power-law), and 10 -11 photons 
cm" 2 s~ 1 for the MAGIC VHE 7-ray observations. Upper limits 
are given at the 95% confidence level, following the prescription by 
IRolke e t al. (2005j). The data were taken at 5.0 GHz for the ra- 
dio observations (except for the MERLIN data in November 2006, 
which were taken at 6.0 GHz), at 0.5—10 keV by Chandra, and at 
0.3-5.0 TeV by MAGIC. 

November included only MAGIC and MERLIN observa- 
tions. 

The range of observed orbital phase in the October 
campaign was O rb = 0.6-0.7, for which the TeV max- 
imum was detected in the past (| Albert et al. I l20"06a|) . 
and or h = 0.44-0.57 in the November campaign. The 
super-orbital phase for both campaigns was </>sup = 0.4. 

2.1. Radio Observations 

We observed LS I +61 303 with several radio interfer- 
ometric arrays during 25-26 October 2006, including the 
Multi-Element Radio Linked Interferometer (MERLIN) 
in the UK, the European VLBI Network (EVN), and the 
Very Long Baseline Array (VLBA) in the USA, all of 
which observed at 5.0 GHz. Additionally, we also moni- 
tored the flux and large angular scale structural and flux 
density variations of LS I +61 303 using MERLIN at 5.0 
GHz on 27 and 28 October 2006, and at 6.0 GHz for the 
period from 16 to 20 November 2006. 

The MERLIN array included six antennas (Defford, 
Cambridge, Knockin, Darnhall, Jodrell Bank (Mkll), 
and Tabley) for most of our observations, yielding syn- 
thesized beams of about 50 to 70 mas (corresponding to 
projected linear sizes of 100 to 140 AU). The EVN ob- 
servations on 26 October 2006, which were the first ever 
carried out for LS I +61 303 using the e-VLBI technique 
(Szomoru 2006), included five antennas spread over Eu- 
rope (Cambridge, Jodrell Bank Mkll, Medicina, Torun, 
and the Westerbork phased array) , yielding a synthesized 
beam of about 7 mas (or 14 projected AU). The data 
were directly streamed to the correlator at JIVE (Joint 
Institute for VLBI in Europe, Dwingeloo, The Nether- 
lands) through the Internet, and this was one of the first 
science observations that achieved 256 Mbps sustained 
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data rate with the e-VLBI technique. The VLBA obser- 
vations on 25 and 26 October 2006 included ten 25-m 
antennas spread over the USA, yielding a beam size of 
4.6x2.1 mas, corresponding to a projected linear reso- 
lution of ~9.2 AU and ~4.2 AU in right ascension and 
declination, respectively. 

All three interferometric arrays observed LS I +61303 
in phase-referenced mode. LS I +61 303 and the 
bright (S5GHz~ 600m Jy), nearby International Celes- 
tial Reference Frame (ICRF) source J0244+6228 were 
alternately observed through each observing run, the 
phases of J0244+6228 being transfered to the position of 
LS I +61 303 in the posterior data analysis. J0244+6228 
also served as amplitude calibrator for our observations. 
We performed standard calibration and data reduction 
within the NRAO Astronomical Imaging Package Sys- 
tem (AIPS; Diamond 1995). We also used standard hy- 
brid mapping techniques within AIPS to obtain the flux 
densities shown in Figure [JJ and Table [TJ and the radio 
images shown in Figures [2] and [3] 

2.2. X-ray Observations 

We obtained Chandra X-ray observations of 
LS I +61303 on 25 October 2006 through the Di- 
rector Discretionary Time program (Chandra Obsld 
8273). The observations were carried out using ACIS-I 
for a total exposure time of 20.0 ks. At the time of 
the observation, LS I +61 303 was expected to be in 
a high-state and its high X-ray brightness could have 
resulted in an excessively high count rate producing 
appreciable pile-up in the observations. In order to 
minimize pile-up effects, LS I +61303 was offset by 
8' from the ACIS-I aim-point, thus smearing its image 
and reducing the count rate at the source peak. We 
further used 1/4 sub-array, reducing the exposure time 
of individual frames from the nominal exposure time of 
3.2 s down to 0.8 s. 

Data reduction of the X-ray observations of 
LS I +61303 have been performed using the Chandra 
X-ray Center software CIAO V3.4. The data reduction 
included the application of standard filters and rejection 
of events with bad grades and those originating from 
bad pixels. The background count rate is consistent 
with the quiescent background (Markevitch 2001), and 
no time intervals of enhanced background needed to 
be removed. The processed observations have a useful 
exposure time of 19.1 ks. Due to the high count rate of 
LS I +61 303 out-of-time events were not negligible and 
produced a noticeable streak along the columns. These 
have been corrected using the CIAO task "acisreadcorr" . 
Finally, notice that as of CALDB v3.1.0 (June 2005), 
the CIAO task "acis_process_events" is routinely used 
by the level 1 processing pipeline to mitigate the charge 
transfer inefficiency (CTI) that affects the Chandra 
ACIS-I front illuminated chips. Therefore, no further 
CTI correction was applied. Lightcurves and spectra of 
LS I +61 303 were obtained using standard CIAO tasks 
and analyzed using HEASARC FTOOLS and XSPEC 
Vll.2.0 routines (Arnaud 1996). 

2.3. VHE ^-ray Observations 

Observations in the very high energy (VHE) 7-ray 
band (E 7 > 100 GeV) were scheduled with the MAGIC 
telescope on 26-28 October and 16-19 November 2006 



(see Table [TJ). These observations are part of an ex- 
tensive, standalone observing campaign c arried out be- 
tween September a nd December 2006 (|Albert et al. I 
12008. in preparation! ) . Bad weather conditions at the 
Observatorio del Roque de los Muchachos, however, pre- 
vented us from obtaining useful data on 26 and 28 Octo- 
ber 2006. 

The observations we re carried out in t he false-source 
track (wobble) mode (|Fomin et al.l Il994) . with two di- 
rections at 24' distance and opposite sides of the source 
direction, which allows for a reliable estimation of the 
background with no need of extra observation time. 

The data were analyzed using t he standard MAGIC 
calibration and a nalysis software (| Albert et al.l l2006bl : 
iGaug et all [2005). Data runs with anomalous event 
rates were discard ed for further analysis. Hillas vari- 
ables (lHillaslll985h were combined into an adimensional 
7/hadron discriminator (hadronness) and an energy esti- 
mator by means of the Random Forest classification algo- 
rithm, which takes into acco unt the correlat i on between 
the d ifferent Hillas variables (jBrciman 2001: Boc k et aTl 
2005]). The incoming direction of the primary 7-ray 
events was estimated using the DI SP method, suited 
for observations with a single IA CT (|Fomin et al.lll994l ; 
iDomingo-Santamaria et ah l2005f ). 

3. RESULTS 

In this section we present the results obtained from 
the different multi- wavelength observations we have per- 
formed, and put them into the context of the past mea- 
surements in the different bands. They are summarized 
in Tableland Figure [TJ 

3.1. Radio Results 

The total radio flux density obtained with MERLIN 
shows a decline between October 26 (</> or ^ = 0.66) and 
28 (0 or K = 0.70), from ~ 35 to ~ 30 mJy, and a peak on 
18 November (</> or b = 0.50) at ~ 80 mJy. At the super- 
orbital phase of the observations (<^sup = 0.4) the radio 
source is in the weak state. The predicted flux of the 
radio flux is between 50-100 mJy at the orbital phases 
^orb ~ 0-9. We therefore measure a flux compatible 
with the predicted one, although at a much earlier phase 
value. However, for the weak state, secondary peaks of 
comparab le flux show up at o ther orbital phase values 
(see, e.g.. lDhawan et aLll2006h . 

In Figure [2] we show the radio images corresponding 
to 25-26 of October, where all three arrays observed si- 
multaneously LS I +61303. In Figure [3] we show the 
MERLIN images obtained on 27 and 28 October 2006 at 
5.0 GHz, and on 16 through 20 November at 6.0 GHz. 
Our MERLIN observations show no evidence for large 
angular scale str uctures, contrary t o what had been pre- 
viously claimed (|Massi et al.ll2004T ). hence excluding the 
existence of persistent jets at these scales. The structure 
of LS I +61 303 genuinely follows the beam shape of the 
MERLIN array at all epochs. Model fits to the u— v plane 
of the MERLIN data show that the projected intrinsic 
size of the radio emitting region of LS I +61303 is, at 
any observing epoch, no larger than ~6 mas (~12 AU), 
as confirmed by our higher resolution images obtained 
with the EVN and the VLBA. 

Our 5 GHz VLBA images on 25 and 26 October show a 
very bright and unresolved component plus an extended 
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Fig. 1. — Radio, X-ray and VHE 7-ray light-curves obtained with VLBA, EVN, MERLIN (left-hand scale), Chandra (right-hand scale 
in 10 — 11 ergs cm -2 s _1 ) and MAGIC (right-hand scale in 10" 11 photons cm -2 s — 1 ) during the two observing periods (October and 
November 2006). The horizontal error bars show the time spanned by the different observations. The shaded are a marks the ran ge of X-ray 
flux values previously reported (see Table [3j . The upper axis shows the orbital phase using the ephemeris from [Gregory ( 2003) . 
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Fig. 2.— Radio images of LS I +61 303 obtained on 25 October 2006 with the VLBA (panel a) and on 26 October 2006 with MERLIN 
(panel b), EVN (panel c), and VLBA (panel d). The synthesized beams of the images (bottom left corner in each panel) are 70x41 
mas (position angle PA=29°), 7.2x5.8 mas (PA=35°), and 4.6x2.1 mas (PA=-81°) for the MERLIN, EVN, and VLBA observations, 
respectively. The projected linear resolution of our VLBA observations is thus of 9.2 AU and 4.2 AU in RA and DEC, respectively. In 
panel a, we also show overlaid the 8 GHz VLBA image on 2 February 2006 (grey scale) convolved with our 5 GHz VLBA beam. This date 
corresponds to the same phase of LS I +61 303 (fli or b ~0.62), and they show a striking similarity. In all cases, the origin of coordinates is 

set at the VLBA peak of brightness on 25 October 2006, and the contours are drawn at (3,3 \/3,9,...) times the off-source rms (see Table 
1). 
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Fig. 3. — Radio images of LS I +61 303 obtained with MERLIN during our October 2006 and November 2006 campaigns. The images 
correspond, from left to right and top to bottom, to observations carried out on 27 and 28 October 2006 at a frequency of 5.0 GHz, and 
on 16, 17, 18, 19, and 20 November 2006 at a frequency of 6.0 GHz, respectively. Radio contours are drawn at (3,3 \/3,9,...) times the 
off-source rms (see Table 1). Note that the shape of the radio brightness distribution of LS I +61303 genuinely follows the synthesized 
beam for each epoch (upper left inset in each panel). 
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radio emission to the east-southeast (see Figure [2]). The 
general as pect of the radi o emiss ion recalls the structure 
found by iDhawan et al71 |2006). In fact, panel (a) of 
Figure [2] shows both our 5 GHz VLB A image on 25 Oc- 
tober 2006 (<^> or b = 0.62, contours) overlaid on top of the 
8.4 GHz VLBA data on 2 February 2006 ((f> orh = 0.62, 
grey scale), kindly provided by V. Dhawan in advance of 
publication. 

The projected distance from the peak of the radio 
brightness distribution to the edge of the extended re- 
gion is, as imaged with the VLBA at 5.0 GHz on 26 
October 2006, of about ~8.2 mas (16.4 AU). If the max- 
imum time to fill in this region with radio (synchrotron) 
emitting particles is less that the time spanned by our 
two consecutive VLBA observations (1 day), the implied 
outflow velocity is at least of 21500 km/s, or w > 0.09 c. 

The coordinates of the peak of radio brightness 
distribution of LS I +61 303 , as obtained from our 
VLBA observations (using task IMFIT), were RA = 

02 h 40 m 31. 6638849 s , DEC = 61°13'45. 592235" on Oc- 
tober 25 and RA = 02 h 40 m 31. 6638756 s , DEC = 
61°13'45. 592496" on October 26, with and estimated ac- 
curacy of 12/xas and 7/xas in RA and DEC, respectively. 
This shift of the brightness peak corresponds to a day- 
to-day projected speed of 904+60 km s^ 1 , in excellent 
ag reement with the typic al value of ~ 1000 km s" 1 found 
bv lDhawan et al.~l (2006) along a complete orbital cycle. 

The peak flux density varied significantly between our 
two VLBA observing runs, decreasing from 24.3 mJy/b 
on October 25 to 15.4 mJy/b on October 26. The total 
5 GHz flux density varied from 39.5+0.2 mJy (radio lu- 
minosity L K = (9.5 ± 0.1) x 10 29 ergs s" 1 ) to 32.7+0.2 
mJy (L R = (7.8 + 0.1) x 10 29 ergs s -1 ). We therefore sug- 
gest that the drop seen in the total radio flux density of 
LS I +61 303 between our two consecutive VLBA obser- 
vations is directly related to the change in the peak flux 
density, which hints towards a physical link with struc- 
tural variations in the innermost ( < 3 mas projected ra- 
dius, or 6 AU at a distance of 2 kpc) region of the source, 
which our observations cannot resolve. 

3.2. X-ray Results 

LS I +61303 is clearly detected in our Chandra ob- 
servation, with a background-subtracted averaged count 
rate of 1.067+0.008 cnts s _1 in the energy band 0.5- 
10.0 keV. This count rate would have indeed resulted in 
considerable pile-up for a standard observational setup. 
The X-ray spectrum (Figure [J) shows a broad peak at 
~1.6 keV and a hard energy tail that extends up to 
9 keV. Assuming a foreground interstellar absorption 
wi th solar abundances and the absorp t ion cro ss-sections 
of iBalucinska-Church fc McCammon I (|1992f) . the X-ray 
spectrum of LS I +61 303 can be reasonably well fit- 
ted by either an absorbed power law, or an absorbed 
Bremsstrahlung model (Fig.rjJ. The best fit parameters, 
goodness of the fit (reduced x 2 ), and implied source flux 
and luminosity in the energy range of 0.5-10.0 keV for 
these two models are listed in Table [2] In both cases, 
the unabsorbed X-ray luminosity of LS I +61 303 in the 
energy range of 0.5-10.0 keV is ~ 10 34 ergs s _1 . 

We have also investigated the short-term variability of 
LS I +61 303 X-ray flux and hardness ratio (defined as 
the ratio of the count rate between 2.0-10.0 keV over that 
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Fig. 4.— Chandra ACIS-I spectrum of LS I +61 303 over- 
plotted with the best-fit absorbed power law (red histogram) and 
Bremsstrahlung (blue histogram) models (see Table [3] for details). 
The lower panel shows the relative residuals of the fits (AI) in terms 
of the bin standard deviation (tr) for both the absorbed power law 
(red dots) and Bremsstrahlung (blue dots) models. The errors 
bars in both the spectra and residual plots are 1-cr. The insets 
show the x 2 plots as a function of the power law index, T, and 
Nn (upper inset), and of the plasma temperature, kT, and TVjj 
(lower inset) of the spectral fits to the Chandra ACIS-I spectrum 
of LS I +61 303 using absorbed power law and Bremsstrahlung 
models, respectively. The contours represent 68%, 90%, and 99% 
confidence levels. 



TABLE 2 

X-RAY SPECTRUM BEST-FIT PARAMETERS a 



Model 


Power-law 


Bremsstrahlung 


X 2 /DoF 


413.95/370=1.12 


428.24/370=1.16 


iV H (10 21 cm" 2 ) 


5.5+0.5 


4.8+0.4 


T or kT 


1.53+0.07 


17±| keV 


r0.5-10.0keV 
1 obs 


1.87+0.16 


1.84+0.05 


£0.5-10.0 keV 


1.18+0.11 


1.11+0.04 



a Best fit parameters to the X-ray spectrum for LS I +61 303 ob- 
tained using absorbed power-law and Bremsstrahlung models. The 
flux between 0.5 and 10.0 keV is expressed in 10 -11 ergs cm -2 s _1 
and the luminosity in 10 34 ergs s — 1 and for an assumed distance 
of 2 kpc. 
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Fig. 5. — Temporal evolution of the Chandra ACIS-I hard- 
ness ratio (top) and count rate in the energy band 0.5-10.0 keV 
(bottom). The bin size is 1,000 s. and the error bars arc l-<r. 
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TABLE 3 

X-RAY OBSERVATIONS OF LS I +61 303 IN THE ~1-10 KEV RANGE a 



Obs. ID 


Date 


Exposure 


^oro 


'/'sup 


r 


/ 


AO 


03-02-1994 


18.0 


0.20 


0.61 


1.71 ± 0.08 


0.72 


Al 


09-02-1994 


19.3 


0.43 


0.61 


1.83 ± 0.08 


0.52 


R0 


01-03-1996 


8.9 


0.79 


0.07 


— 


0.84 


Rl 


04-03-1996 


9.3 


0.91 


0.07 


- 


1.21 


R2 


07-03-1996 


9.1 


0.03 


0.07 


- 


0.69 


R3 


10-03-1996 


9.1 


0.11 


0.07 


- 


0.65 


R4 


13-03-1996 


10.0 


0.22 


0.07 


- 


0.96 


R5 


16-03-1996 


9.5 


0.33 


0.07 


- 


1.32 


R6 


18-03-1996 


9.5 


0.42 


0.08 


— 


2.00 


R7 


24-03-1996 


11.0 


0.64 


0.08 


— 


1.22 


R8 


26-03-1996 


13.7 


0.71 


0.08 


- 


1.03 


R9 


30-03-1996 


14.6 


0.87 


0.08 




0.92 


SO 


22-09-1 997 


12.0 


0.31 


0.41 


J-.VJO _l_ U.1U 




SI 


26-09-1997 


8.6 


0.45 


0.41 


1.56 ± 0.09 


1.40 ±0.01 


XO 


05-02-2002 


6.4 


0.56 


0.37 


1.60 ±0.02 


1.39 ±0.01 


XI 


10-02-2002 


6.4 


0.76 


0.37 


1.53 ±0.02 


1.35 ±0.01 


X2 


17-02-2002 


6.4 


0.01 


0.37 


1.74 ±0.05 


0.50 ±0.01 


X3 


21-02-2002 


7.5 


0.18 


0.38 


1.57 ±0.02 


1.36 ±0.01 


XI 


16-09-2002 


6.4 


0.97 


0.50 


1.60 ±0.07 


1.39 ±0.01 


X5 


27-01-2005 


48.7 


0.60 


0.02 


1.62 ± 0.01 


1.29 ±0.01 












1.83 ± 0.01 


0.40 ±0.01 


Co 


07-04-2006 


49.9 


0.03 


0.28 


1.25 ± 0.09 


0.71 ± 0.18 


Cl 


25-10-2006 


20.0 


0.61 


0.40 


1.53 ± 0.07 


1.87 ± 0.02 


a The first 


column shows an ID 


composed of a letter 


standing for 


the X-ray satellite (A 


= ASCA ILeahv et al. 


1997), R = 



RXTE JHarrison et al. I [20001) S = BeppoSAX (ISidoli et al.ll2006T) , X = XMM-Newton flSidoli et al.l 12001 IChernvakova et al .112001) . C 
= Chandra (Parcdcs ct al. 2007, and this work)) and an ordering index. Observation Cl corresponds to the one performed during our mul- 
tiwavelength campaign. The exposure is in ks. (frorb ant ^ 'Psup are the orbital and super-orbital phases of the beginning of the observation, 
computed using To = 43366.275, P or b = 26.4960 days and Psup = 1667 days (Gregory 2002). The photon index (r) is obtained from a 
fit to an absorbed power-law model. The flux (/) is expressed in 10 -11 ergs cm -2 s _1 , and is integrated between 2-10 keV for RXTE, 
BeppoSAX and XMM-Newton, between 0.3-10 keV for observation CO and between 0.5-10 keV for Cl. 



between 0.5-2.0 keV), shown in Figure [5] The Chandra 
ACIS-I count rate of LS I +61 303 in the energy band of 
0.5-10.0 keV varied by a 25%, from - 1.25±0.03 cnts s" 1 
down to ~ 1.00 ± 0.03 cnts s^ 1 , within a timescale of 
1-2 hours, showing a clear decline in the final third of 
the observation. The probability for this variation to be 
a statistical fluctuation of a constant flux is less than 
10~ 15 , as obtained from a x 2 fit of a constant function to 
the data set. Fast flux variations of up to a 70% within 
few ho urs time-scale have been reported by ISidoli et"ai1 
(2006). In their case, these were accompanied by varia- 
tions in the hardness ratio (HR = 2-12 keV/0.3-2 keV) 
of about 30%, which our data do not confirm. It must 
be noted, however, that assuming a linear correlation be- 
tween flux and HR, we should observe less than ~ 10% 
variations in the HR during our observation, which is 
below our statistical error. 

Table [3] summarizes the results of the observations of 
LS I +61303 in the ~1-10 keV energy range existing 
in the literature, including that obtained in our Chandra 
observations. We measure a flux of 1.9 x 10~ n ergs cm~ 2 
s _1 and a photon index 1.53. The historical values ob- 
tained over several years range between [0.4 — 2.0] x 10~ n 
ergs cm~ 2 s _1 for the flux and [1.25 — 1.83] for the photon 
index. Therefore, we observed the source in a particu- 
larly high and hard state. This is shown also in Fig- 
ure El where we plot the historical values of the photon 
index against the flux, in the soft X-ray range. The Pear- 
son's correlation coefficient between these two quantities 
is r = —0.46, but becomes r = — 0.91 when the outlier 
point obtained bv lParedes et al.l (|2007f ) is removed from 
the computation. This is, to our knowledge, the first time 
that this correlation is shown, and confirms the claim by 
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Fig. 6. — Correlation between the X-ray flux and photon index 
for the existing observations of LS I +61 303 in the ~1— 10 keV 
energy band, obtained by using public data from the literature 
and our own data. (See main text for details). 

ISidoli et al.l (|2006l ) that the source is harder when it is 
brighter. 

3.3. VHE j-ray Results 

LS I +61 303 was detected with MAGIC only on Octo- 
ber 27 (^oj-v, = 0.66), with a significance of 4.5cr. For the 
rest of the nights, no detection above 2cr was found, and 
we derived the corresponding upper limits to the integral 
flux (see Table P. 

On October 27, the measured average flux above 300 
GeV corresponds to 15% of the Crab Nebula flux at 
these energies. The VHE 7-ray source is point-like for 
the MAGIC angular resolution (0.1°) and the location is 
compatible with that of LS I +61 303 . The energy spec- 
trum is well fitted by an unbroken power-law with index 
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Fig. 7. — Quasi-simultaneous LS I +61330 spectrum including 
radio (VLBA, open circles), X-rays (Chandra, squares) and VHE 
7-rays (MAGIC, stars) data from the period 25-26 October 2006, 
along with the data from November 2006: average MERLIN flux 
density (filled circle) and upper limits from MAGIC (arrows). 

a = — 2. 7± 0.5 ±0.2, where the quoted errors correspond 
to the statistical and systematic uncertainties, respec- 
tively. No significant variations of the absolute flux were 
detected within that night. 

Previous observations of this source with the MAGIC 
(lAlbert et al. I l2006ai 12008. in preparation! ) and VERI- 
TAS (|Acciari et al. II2008D telescopes have shown a point- 
like source, with a peak of ~ 15% Crab flux inten- 
sity at orbital phase 4>orb ~ 0.65, and a spectral index 
a ~ —2.6, hence in agreement with the results derived 
from our observations. 

4. DISCUSSION 

The comparison between VHE 7-ray and radio data 
during the October campaign (see Figure [1} shows a de- 
tection at TeV energies for or ^ = 0.66 with a flux level 
of ~ 15% of the Crab Nebula flux, during a period when 
the radio emission is constant at 35 mJy. During the 
November campaign, however, there is no detection at 
TeV energies, while the r adio data show a peak flux twice 
as high as in October. lAlbert et al. I (|2006a| ) reported 
on radio and TeV peaks detected almost simultaneously, 
while for our campaing we see the TeV peak for a flat 
and low radio flux and a radio peak for no significant 
TeV emission. Therefore, we exclude a general TeV- 
radio correlation. A plausible explanation is that the 
emissions are produced by different particle populations. 
On the other hand, the detections at X-ray and TeV en- 
ergies, both at particularly high flux values and one day 
apart, might point to a correlation between X-ray and 
TeV fluxes, and hence to the fact that both radiations 
are produced by the same population of particles. Yet, 
our data are too scarce to make any solid conclusion at 
present. 

The results obtained by radio imaging at different an- 
gular scales show that the size of the radio emitting re- 
gion of LS I +61 303 is constrained below ~6 mas (~12 
projected AU), and the presence of persistent jets above 
this scale is therefore excluded. 

We resolve a radio-emitting region at 5.0 GHz with 
VLBA, extending east-southeast from the brighter, un- 
resolved emitting core. The outflow velocity implied by 
our observations is at least of ~0.1 c, in agreement t o 
what previously suggested bv lChernvakova et al.l (|2006|) . 



for an over-pressured puls ar wind, although rece nt hy- 
drodynamical simulations (Bogovalo v et al.ll2008h show 
that the shocked pulsar wind could be relativistic. 

The comparison of our VLBA image from 25 Octo- 
ber with an image at 8.4 GHz obtained at a similar 
orbital phase but 10 orbital cycles apart shows a high 
level of similarity bet ween them. It is worth noting that 
iDhawan et all 0006) obtained high-resolution radio im- 
ages of the source along a complete orbital cycle, and 
found an extended feature whose orientation with respect 
to the central core varied by 360° along the orbit. There- 
fore, the similarity in the morphology and orientation 
between both images suggests a significant level of pe- 
riodicity and stability of the physical processes involved 
in the radio emission. This result points to the fact that 
the extended radio emission is produced by the inter- 
action of steady flows (from a relativistic pulsar wind, 
jet and/or stellar wind) rather than by the interaction 
of such an outflow with wind clumps. We note that if 
the radio emission is produced by a milliarcsecond scale 
jet, the required stability and periodic behavior of such a 
jet is difficult to reconcile with the non-persistent nature 
of a large scale (~100 mas), putative relativistic jet, as 
deduced from our ME RLIN observa t ions i n combination 
with those obtained by iMassi et alj (|2004[) . 

Finally, we combine our X-ray data from 25 October, 
the TeV data from 26 October and the average VLBA 
in both nights to produce a quasi-simultaneous multi- 
wavelength spectrum, including radio, X-ray and VHE 7- 
ray observations, which we display in Figure[71 However, 
the simultaneous data cover only a small range of the or- 
bital phase of LS I +61 303 . Given the high variability 
of the physical conditions of the system along the or- 
bit, more simultaneous, multiwavelength data, and par- 
ticularly involving longer exposure times, orbital phase 
coverage and redundancy, will shed further light in our 
understanding of this peculiar object. 
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